Draft version August 25, 2011 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



A WISE VIEW OF STAR FORMATION IN LOCAL GALAXY CLUSTERS 

Sun Mi Chung 1 , Peter R. Eisenhardt 2 , Anthony H. Gonzalez 1 , Spencer A. Stanford'*, Mark Brodwin 4 , Daniel 

Stern 2 , Thomas Jarrett 5 

Draft version August 25, 2011 



o 

(N 
< 

m 
(N 



ABSTRACT 

We present results from a systematic study of star formation in local galaxy clusters using 22/im 
data from the Wide- field Infrared Survey Explorer (WISE). The 69 systems in our sample are drawn 
from the Cluster Infall Regions Survey (CIRS), and all have robust mass determinations. The all-sky 
WISE data enables us to quantify the amount of star formation, as traced by 22/im, as a function 
of radius well beyond R200, and investigate the dependence of total star formation rate upon cluster 
mass. We find that the fraction of star-forming galaxies increases with cluster radius but remains 
below the field value even at 3R2oo- We also find that there is no strong correlation between the 
mass-normalized total specific star formation rate and cluster mass, indicating that the mass of the 
host cluster does not strongly influence the total star formation rate of cluster members. 
Subject headings: galaxies: clusters: general — galaxies: evolution 
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1. INTRODUCTION 

It is well-established that the fraction of star-forming 
galaxies declines as a function of increasing local galaxy 
density in the low redshift universe . Also known as 
the star formation -density relation (|Lewis et al.1 120021 : 
iGomez et ah! 120031: iBalogh et all |200j), this correlation 
has been confirmed in many studies, primarily using op- 
tical and UV data to trace star formation in massive 
galaxy clusters and field environments. Mid-infrared 
data from the Infrared Satellite Observatory (ISO) and 
the Multi- band Imaging Photometer for Spitzer (MIPS) 
have also revealed the presence of highly obscured, dusty 
star forming galaxies, previously undetec t ed by opti- 
cal o r UV surveys (e.g. ICedres et al.l 120091: ISmith et al.l 
12010ft . While the sensitivity of MIPS has enabled de- 
tailed studies of obscured star form ation in individua l 
local and distant galaxy clusters (e.g . IGeach et alll200 
Marcillac et al.l 120071: IDressler etaLll2009t [Haines et al 



20091: iChung et all 120101 IFinn et all 12010) . there are 
still only a small number of low redshift clusters that 
have been systematically surveyed for dusty star-forming 
galaxies out to the virial radius. 

There remain many uncertainties in the relationship 
between star formation in clusters and their global clus- 
ter properties. In particular, several studies have tried 
to understand the correlation between cluster mass and 
the mass-normalize d cluster s t ar for mation rate (SFR). 
While results from iBai et all (J2009) suggest that there 
is no strong correlation between c luster speci fi c SFR 
and cluster mass, ot hers such as IFinn et al.l 12005). 
iPoggianti etaTJ ([2006), and lKovama et al.l (|2010) argue 
that cluster specific SFR decreases with cluster mass. 
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The large spatial coverage required to observe dusty 
star-forming galaxies in low redshift clusters out to the 
cluster infall regions has thus far hindered our ability 
to understand how star formation is affected by global 
cluster properties such as cluster mass. In this paper 
we exploit dat a from the Wide-fie ld Infrared Survey Ex- 
plorer (WISE; IWright et all 12010) to overcome this ob- 
servational challenge and present results on obscured star 
formation and how it relates to cluster mass and radius 
out to 3R200 in a sample of 69 clusters at z < 0.1. R200 
and M200 are commonly used interchangeably with virial 
radius and total cluster mass, respectively. R200 is the 
radius within which the average density is 200 times the 
critical density of the universe and M200 is the mass en- 
closed within that radius. 

2. DATA 
2.1. WISE 

WISE is a medium-class explorer mission funded by 
NASA and has completed observations of the entire sky 
in four infrared bands: 3.4, 4.6, 12, and 22/im (Wl to 
W4, respectively). WISE scanned the sky with 8.8 sec- 
ond exposures, each with a 47 arcmin field of view, pro- 
viding at least eight exposures per position on the ecliptic 
and increasing depth towards the ecliptic poles. The in- 
dividual frames were combined into coadded images with 
a pixel scale of 1.375 arcsec per pixel. Cosmic-rays and 
other transient features were removed via outlier pixel 
rejection. 

The photometry used for our analyses is point spread 
function (PSF) fitted magnitudes from the "first-pass op- 
erations coadd source working database" created by the 
WISE data reduction pipeline. Galaxies in our cluster 
sample have a diffraction limited resolution of 12" (full 
width half maximum) in the 22/xm band. We have con- 
firmed from W4 coadded images that all star-forming 
galaxies considered in our analyses appear unresolved in 
the 22/zm band, and have PSF photometry reduced % 2 
values less than 1.5. Therefore we use the PSF magni- 
tudes from the first-pass photometric catalog to obtain 
estimates of total flux. 
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For the minimum coverage of 8 overlapping frames, the 
sensitivity for S/N = 5 in the W4 band is 6 mJy, includ- 
ing u ncertainty due to source confusion (Wrig ht~et al.l 
2010). To ensure an unbiased comparison of global SFRs 
and total IR luminosities of clusters at different redshifts, 
we impose a lower limit of SFR=4.6 M Q yr _1 on our en- 
tire cluster sample, which is equivalent to a total IR lu- 
minosity of Lir> 4.7 x 1O 1O L , and corresponds to the 6 
mJy flux limit at z = 0.1. We hereafter refer to our sam- 
ple of star-forming galaxies as demi-LIRGs, which have 
nearly half the total IR luminosity of a luminous infrared 
galaxy or LIRG. However we note for future extragalactic 
studies using WISE data that most coadded observations 
will have at least 12 overlapping frames and hence better 
sensitivity than the conservative 6 mJy limit we adopt 
in this paper. Additional information regarding WISE 
data processing is available from the Preliminary Data 
Release Explanatory Supplement. Q 

2.2. Cluster Sample 

We use the Cluster Infall Regions (CIRS; Rines & Di- 
aferio 2006) sample because it provides high-fidelity mass 
estimates, is at sufficiently low redshift to enable detec- 
tion with WISE of strongly star-forming galaxies, and 
has extensive spectroscopy for membership determina- 
tion. The CIRS sample consists of 72 low-redshift X-ray 
galaxy clusters identified from the ROSAT All-Sky Sur- 
vey that are within the spectroscopic footprint of SDSS 
Data Release 4. The redshift range of the CIRS clusters 
is 0.003 < z < 0.1, with a m edian of z = 0.06 . Clus - 
ter masses are available from IRines fc Diaferid {2006), 
who utilize the caustics infall patte rn to determine tota l 
dynamical cluster mass and R 200 piaferio et alj|2005l) . 
The clusters in this paper consist of the entire CIRS 
sample, excluding three clusters at z < 0.006, which 
leaves 69 remaining clusters with a minimum redshift of 
z = 0.0204. 

2.3. SDSS DR7 

Optical photometric and spectroscopic data are ob- 
tained from the Slo an Digital Sky Survey Data Re- 
lease 7 (SDSS DR7) (|Abazaiian et alj I2009ji . which are 
90% spectroscopically complete for galaxies with r < 
17.77 and half-light surface brightness ^50 < 24.5 mag 
arcsec -2 . However, the spectroscopic completeness is 
lower in high-density regions such as in the core of galaxy 
clusters, due to fiber collisions. Adjacent fibers cannot 
be placed closer than 55 arcsec from each other, which 
corresponds to a separation of 63 kpc at z = 0.06. To 
verify that the spectroscopic incompleteness arising from 
constraints on fiber placement has a negligible effect on 
radial trends of star formation in clusters, we look at the 
fraction of W4-bright sources as a function of distance 
from the cluster center. 

Among the r < 17.77 sample, the spectroscopic com- 
pleteness of W4-bright sources (W4 S/N > 5) is —70% 
within the central 5 arcmin radial bin, and increases to 
—80% in the 5-10 arcmin bin. In other words, we are 
not missing a significant fraction of star-forming galax- 
ies in the cluster core due to fiber collisions. The spec- 
troscopic completeness for all galaxies (regardless of IR 

http://wise2.ipac.caltech.edu/docs/release/prelim/index.html 



emission) with r < 17.77 shows a similar behavior with 
distance from the cluster center, and on average reaches 
at least —80% completeness beyond the central 5 arcmin 
radius. Since the spectroscopic incompleteness of both 
W4-bright and optically bright galaxies reaches —80% 
beyond the cluster core, our results on radial trends of 
star formation are negligibly affected by spectroscopic in- 
completeness. In addition, we are measuring trends on 
scales of hundreds of kpc, and therefore not significantly 
affected by small differences in spectroscopic complete- 
ness near the cluster core. 

We use photometric data from SDSS DR7 to estimate 
stellar masses for cluster members, using the tight cor- 
relation between stellar mass-t o- light (M/L) ra tio and 
optical colors, as determined bv lBell et al.l (|2003l ). With 
a sample of more than 10,000 opti cally bright galaxies 
(13 < r < 17.5), IBell et al.l (|200l construct a grid of 
stellar population models with a range of metallicities 
and star formation histories, then compare the best-fit 
galaxy templates with evolved zero redshift templates to 
determine present-day M/L ratios. We use the relation 
between g—r color and r-band stellar M/L ratio to derive 
stellar masses for our sample. The estimated uncertainty 
of the color-based stellar M/L r atios, including random 
and systematic errors, is -45% (pellet al.M2003h . 

3. ANALYSIS 

3.1. Cluster Membership 

Spectroscopic redshifts from SDSS DR7 are used to 
determine membership for each of our clusters. Figure [1] 
illustrates an example of using the caustics infall method 
to determine cluster membership for Abell 1377, the most 
massive cluster in our sample. It shows the difference in 
radial velocity (cz) of galaxies in the cluster field with 
respect to the cluster systemic velocity, as a function 
of projected distance from the cluster center. Galaxies 
that are dynamically bound to the cluster form a well- 
defined region that decreases in velocity offset as a func- 
tion of projected radius. We define the galaxies that are 
within the edge of this envelope and within a projected 
radius of <3R2oo, as cluster members. While the aver- 
age turnaround radius for the CIRS sample is — 5R200, 
we restrict our cluster galaxy sample to within 3R200, 
since cluster infall patterns are generally better defined 
closer in to the cluster center. 

Only galaxies with SDSS spectroscopic redshifts are 
considered in the following analyses. We also limit all 
cluster and field galaxies to be brighter than M r = —20.3, 
which corresponds to the 90% spectroscopic complete- 
ness limit of r = 17.77 at z = 0.1. 

3.2. Star Formation Rates and Infrared Luminosities 

Star formation rates and total infrared luminosities 
(Lir) are determined from the WI SE 22/xm photome - 
try, using the relations presented in Rickc et al.l (120091) 
that a re calibrated for MIPS 24^tm data. iRieke et al.l 
( 2009) constructed model average spectral energy distri- 
bution (SED) templates from a sample of local LIRGs 
and ULIRGs and derive correlations between 24//m flux, 
SFR, and Lir . The flux-SFR relation given in their equa- 
tion 1 4 can be used for b oth MIPS 24/zm and WISE 22^m 
data (jRieke et al.ll2009f) . 

The similarity o f the two bandpasses are confirmed by 
I Goto et"ail (|2011[ ) who quantify the correlation between 
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Fig. 1. — Velocity offset versus projected distance from the cluster 
center of Abell 1377, for all SDSS galaxies brighter than r < 17.77 
with spectroscopic redshifts, within a projected 5 Mpc radius of the 
cluster center. Filled circles indicate the cluster members within 
3R.200- Cross symbols (blue) highlight field galaxies that are at 
least 5000 km s — 1 away from the cluster systemic velocity. Red 
and green star symbols represent demi-LIRGs in the cluster and 
field populations, respectively. Open circles are galaxies that are 
included in the SDSS spectroscopic sample but are not chosen as 
cluster members nor field galaxies. 

Lir and the inferred MIPS 24/zm and WISE 22/xm lu- 
minosities. The MIPS and WISE luminosities are in- 
ferred by applying a color-correction to the 18/im AKARI 
flux, an d Lir is derived from fitting the IR SED tem- 
plates of iCharv fc Elbad (|200lD to all six AKARI band- 
passes (9, 18, 65, 90, 140, and 160 /im) for ^600 galax- 
ies at z < 0.1. The resulting correlations between the 
Lir — Lmips24 and Lir — Lwise22 are nearly identical, 
with only a ~4% offset. Therefore we proc eed by using 
the flu x-SFR and Lir — LmipS24 relations of iRieke et al.1 
(2009), assuming the calibration determined from MIPS 
24fim data. 

The total error attributed to the SFR is ~0.2 dex, 
which is dominated by scatter in the -Lm.tpS24-Lir re- 
lation, and does not include uncertainties inherent in the 
assumed st ellar initial mass f unction (IMF). T he IMF 
adopt ed by IRieke et al.1 (120091) is similar to the iKroupal 
(|2002ft and lChabrierl (12603) IMFs, which have relatively 
fewer low mass stars compared to a Salpeter IMF, and 
is more applicable for extragalactic star-formi ng regions 
(|Rieke et alJl!993t lATonso-Herrero et al.ll2001| ). 

3.3. Exclusion of AGN 

A robust estimate of the cluster star formation rates re- 
quires that we first identify and exclude all 22[im sources 
that are dominated by AGN rath e r than star formation. 
As demonstrated in iChung et al.l (|2010D . even a few IR- 
bright AGN can significantly bias the inferred global star 
formation rate of a cluster. There are several methods 
by which we can identify AGN in our data set. The 
first method relies upon WISE colo rs. Similar to th e 
AGN wedge in Spitzer/IRAC data (jStern et al.ll2005h . 
AGN-dominated sources will have red colors in W1-W2 
([3.4]- [4.6]). Specifically, we use the criteria Wl-W2>0.5 
(Vega) to identify candidate AGN in the WISE data 
set. This color selection is similar to, though slightly 



bluer than, the AGN selection determined in Jarrett et 
al. 2011 (submitted) and Stern et al. 2011 (in prep). 
While this single color cut is in general less robust than 
the full AGN wedge, at the low redshifts that are the 
focus of this work, star-forming galaxies should have col- 
ors uniformly blueward of this threshold. We illustrate 
in Figure [5] that our Wl-W2>0.5 criterion works well in 
selecting out AGN from star-forming galaxies. Profile- 
fit photometry is used to determine the W1-W2 colors 
and AGN exclusion. While some of the nearest galaxies 
will be resolved in these bands, we are interested only 
in the W1-W2 color, rather than in single-band photom- 
etry. Comparison of AGN selection based on W1-W2 
color from aperture photometry and profile-fit photome- 
try show that the two methods are similarly successful in 
isolating AGN. However, the color cut based on profile-fit 
magnitudes has a slightly better overlap with the opti- 
cally detected AGN, discussed below. 

A second approach to AGN identification is use of op- 
tical spectroscopy to identify sources that lie in the AGN 
region of the Baldwin-Phill ips- Terlevich (BPT) diagram 
(BPT; IBaldwin et aL|[l98l . Here we use the BPT dia- 
gram as a cross-check on the WISE color selection for the 
subset of sources. Figure [5] shows the emission line ratios 
of [OIIIJ/H^ and [SII]/H Q , for a sample of 136 cluster 
members with Lir > 4.7 x 10 10 L^, of whi c h 27 have 
Wl-W2>0.5. Boundaries from iKewlev et al.l (|2006l) are 
shown as dotted lines and separate regions where nar- 
row emission lines arise from the presence of Seyfcrts, 
LINERs (low ionization narrow emission line regions), 
and HII regions (star-forming galaxies). Galaxies that 
are selected as AGN candidates based purely on having 
a red WISE color (Wl-W2>0.5) are indicated as filled 
(red) circles. 

Figure [5] shows that out of 22 optically identified 
Seyferts, fourteen (65%) are also identified as AGN based 
on the WISE Wl-W2>0.5 selection. Sources flagged as 
quasars (or QSOs) in the SDSS catalog are highlighted 
with (purple) triangle symbols, with all ten quasars being 
independently identified as AGN using the WISE color 
selection. Since the BPT diagnostic is meant to classify 
galaxies based only on narrow emission line ratios, it is 
not surprising that nearly all of the SDSS quasars are not 
properly diagnosed in the BPT diagram. The Ha and 
H/3 emission lines in eight out of the ten SDSS quasars 
have broad wings relative to their neighboring [OIII] or 
[SII] emission lines. Overall, we conclude from Figure [2] 
that the WISE color selection is an effective method of 
excluding AGN from our sample of bright W4 sources, 
identifying 85% of the optically detected Seyferts and 
quasars. 

4. RESULTS & DISCUSSION 

The two central questions that we aim to address are 
how the mean specific star formation rate of cluster 
galaxies (mSSFR) depends upon location within a clus- 
ter, and how the total integrated star formation rate of 
a cluster depends upon cluster mass. Specifically, the 
mSSFR is defined as the total star formation rate in clus- 
ter galaxies at a given projected radius as inferred from 
22/im photometry, divided by the total stellar mass of all 
cluster galaxies at that radius. 

For probing the total integrated star formation rate we 
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Fig. 2. — The BPT diagram for 136 cluster galaxies with 
L IR > 4.7 x 10 10 L w . Dotted curves indicate boundaries from 
Kcwlcv ct al. (2006) that separate galaxies whose emission lines 
originate from Seyferts, LINERs, and HII regions. Solid circles 
(red) highlight galaxies that are identified as AGN based on having 
a red WISE color (Wl-W2>0.5), and triangles (purple) represent 
sources flagged as quasars in the SDSS catalog. 



also define the integrated cluster quantity 
SFR{r < 3i? 20 o) 



cSSFR = 



2nd 



(1) 



which is a useful mass-normalized measure of the total 
star formation rate within the infall region, and referred 
to as the cluster specific SFR (cSSFR). 

Among the 69 clusters, a total of 136 demi-LIRGs are 
detected within 3R2001 of which 27 are determined to 
be AGN based on their W1-W2 color. In the follow- 
ing sections, all SFR quantities are determined from the 
remaining 109 star-forming demi-LIRGs. Eight Seyferts 
identified from the BPT diagram are included in the clus- 
ter star-forming galaxy sample because their W1-W2 col- 
ors are not indicative of AGN activity, and we prefer to 
maintain a uniform WISE selection of AGN among all 
field and cluster galaxies. We note that all results pre- 
sented in this paper are negligibly affected by the exclu- 
sion/inclusion of these eight Seyferts. The coordinates, 
rcdshift, W4 magnitude, SFR, Lir, projected distance 
from cluster center, and AGN flags are listed for the 136 
demi-LIRGs in Table 1. 

4.1. Radial Dependence of Star Formation 

There is a long history in the literature demonstrat- 
ing the existence of a strong radial depend ence for star 
formation (or color) in galaxy clusters (e.g. iLewis et al.1 
2002: IG6mez!t^l27MlBalogh et alJl200l . The main 
strengths of the current analysis are the uniform 22/im 
WISE data, spectroscopic completeness, and existence of 
R200 measurements for the full sample, which together 
provide us with an infrared-based view of star forma- 
tion for a homogeneous sample extending well beyond 
the virial radius. A common approach in the literature 
has been to look at the radial dependence of the frac- 
tion of star-forming galaxies. Here we investigate both 
the dependence of the star-forming fraction, and also the 
mSSFR out to 3R 2 o . 
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Fig. 3. — Mean specific SFR (mSSFR) of spectroscopically con- 
firmed cluster members as a function of projected radius, with 
numbers indicating how many star forming galaxies with Lm > 
4.7 X 1O 1O L0 are in each bin. Dashed line indicates the mSSFR of 
the field sample. 

For comparison, the identical quantities for a field pop- 
ulation are calculated, with the field sample chosen from 
a catalog of galaxies located within a projected radius 
of 5 Mpc from the cluster center, at the cluster redshift. 
Among these galaxies, we choose those with a radial ve- 
locity greater than 5000 km s _1 away from the cluster 
systemic velocity, redshift z < 0.1, and absolute r magni- 
tude M r < —20.3. These are the same redshift and mag- 
nitude limits of the cluster galaxy sample. The highest 
velocity dispersion of our cluster sample is ~960 km s _1 , 
which means that the chosen field galaxies are more than 
5<7 away from the cluster redshift. Figure [T] shows delta- 
velocity (cz) versus projected distance from the cluster 
center for galaxies within a 5 Mpc projected radius of 
Abell 119. The field galaxies indicated with (blue) cross 
symbols are clearly not associated with the cluster galax- 
ies, which appear distinctly confined to a trumpet-shaped 
region. By gathering field galaxies within a projected 5 
Mpc radius from 69 different regions of the sky, we have 
compiled a large enough sample to obtain a representa- 
tive field value of mSSFR and the fraction of star-forming 
galaxies with Lib, > 4.7 x 10 10 L e . There are a total of 
11180 field galaxies, of which 566 are demi-LIRGs. 

Figure [3] shows the mSSFR as a function of r/R2oo 
for our ensemble of 69 clusters, including star formation 
only for galaxies with L IR > 4.7 x 10 10 L (SFR> 4.6 
M Q yr _1 ). We emphasize that by construction our SFR 
limit therefore means that the observed mSSFR is a lower 
bound - but a consistent lower bound across the sample. 

As expected, the mean specific SFR increases with pro- 
jected radius, with the most central bin containing the 
fewest star-forming demi-LIRGs - twelve out of a total 
of 109 demi-LIRGs found in 69 clusters. The mSSFR 
displays a steep increase from the central bin to R200! 
then continues to increase monotonically and nearly flat- 
tens out below the field value at larger radii. This is 
consistent with a low redshift study of Ha star-forming 
galaxies by ILewis et"al"1 (|2002f) , who found that the me- 
dian SFR (normalized by L*) of cluster galaxies reaches 
the field value beyond 3R2oo- 
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Fig. 4. — Ratio of star-forming galaxies to all spectroscopic clus- 
ter members as a function of projected cluster radius, with the 
typical field value shown as a dashed line. 

Such a radial trend can be driven by two factors - 
an increase in the SSFR of the sub-population of star- 
forming galaxies, or an increase in the fraction of star- 
forming galaxies with radius. We find that there is no 
statistically significant change in the SSFR of individ- 
ual star-forming galaxies with radius, implying that the 
trend is driven primarily by a radial gradient in the star 
forming fraction. 

This can be seen in Figure 21 where we directly plot the 
fraction of star-forming galaxies as a function of radius. 
More than one third of the star-forming galaxies reside 
at 2R2oo<r<3R2oo- What is striking is that even at these 
large radii the mSSFR and star-forming fraction remain 
below the field values. In Figure HI the demi-LIRG frac- 
tion is higher in the field relative to the cluster population 
at 3R200 by a factor of ~1.5. One possible interpreta- 
tion of this result is that these are infalling galaxies that 
may already be "pre-processed" in intermediate-density 
environments, such as galaxy groups or filaments (e.g. 
Zabludoff & Mulchaey 1998), which would explain why 
they are suppressed relative to a true field population. 

There is both observational and theoretical support 
in favor of th e suppressed s tar-f orming galaxy fraction 
out to 3R2oo- iBalogh et al.1 (|1998| ) found that the mean 
cluster SFR, which they derived from [Oil] emission of 
galaxies in 15 clusters, is lower than the mean field SFR 
by more than a f acto r of two at 2R 2 qq. Results from 
iLewis et all (|2002t ) and lvon der Linden et ail (|2010D sim- 
ilarly conclude that cluster star formation is suppressed 
at several times the virial radius relative to the field. 

In addition to galaxies being pre-processed at large 
radii, another explanation for the suppressed demi-LIRG 
fraction in Figure 2] could be the contribution of a 
"backsplash" popul ation beyond R^nn- Sim ulations from 
IBalogh et all pOOOT ) and iGill et all^OOl show that up 
to ^50% of galaxies currently located between R 2 no and 
2R200 may have previously traveled inward of the virial 
radius of the cluster. Therefore the suppressed demi- 
LIRG fraction at large radii may, to some degree, reflect 
quenched star formation that occurred within the cluster 
core. However our galaxy sample, which has an abso- 



lute magnitude limit of M r < —20 .3, is sensitive t o mas - 
sive galaxies brighter than M* + 1 iBlanton et al.l (|2003f) . 
whereas backsplash galaxies are on average significantly 
less m assive than first infall galaxies in the cluster out- 
skirts (jGill et al.l 1200 5) . As such, our sample is biased 
against the lower mass backsplash galaxies relative to 
first infall galaxies. Nonetheless, a contribution of back- 
splash galaxies may indeed be partially responsible for 
the suppressed star formation activity at large distances 
from the cluster core. 

4.2. Cluster Specific SFR versus Cluster Mass 

Total mass is, both for astrophysics and cosmology, 
the most fundamental physical galaxy cluster parameter. 
There are multiple physical processes that depend on the 
depth of the cluster potential well (e.g. ram pressure, ha- 
rassment, tidal interactions) which can significantly alter 
the morphologies, gas c ontent, and star formation rates 
in cluster galaxies fe.g. iGavazzi et al.l 2001tjBekki et al.l 
23 
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[20021: lOwen et all 120051: lHavnes et al.l I2007T T Further- 
more, it has been observed that while the total baryon 
fraction (stellar and gas content) remains roughly con- 
stant with cluster mass, the stellar mass decreases and 
gas content in creases in galaxy clusters as a function of 
cluster mass JlGonzalez et a J 120071: iGiodini et all l200l 
lAndreor] I2010T K This implies that the integrated star 
formation efficiency of the history of a cluster is directly 
tied to cluster mass. 

In this section we examine the relation between cluster 
mass and total current cluster star formation to better 
understand how strongly the present-day star formation 
rate depends upon cluster mass. Our approach is to com- 
pute the cluster SFR normalized by cluster mass (cSSFR; 
equation 1) as a function of cluster mass. For each clus- 
ter, the cluster SFR is the sum of SFRs for all member 
galaxies that have SFR>4.6 M Q yr _1 and lie within a 
projected radius of 3R2oo- 

The left panel of Figure shows the cSSFR for 62 
clusters, which consists of the sample of 69 clusters, ex- 
cluding seven clusters that have incomplete spatial cover- 
age in the SDSS DR7 spectroscopic survey within 3R200 • 
The dotted curve shows the limiting detectable cSSFR 
as a function of M200 corresponding to a single demi- 
LIRG with a SFR of 4.6 Mgyr" 1 in a cluster. Clusters 
that have no members above the SFR>4.6 M^yr -1 limit 
are assigned a cSSFR of zero. The binned data (large 
squares) indicate the average cSSFR in each M200 bin, 
including clusters with a cSSFR of zero. Thus the cSSFR 
of the lowest M200 bin is strongly biased by incomplete- 
ness, since nearly all of the clusters at low mass fall be- 
low the SFR>4.6 M yr^ 1 limit and have a cluster SFR 
set to zero. To ensure that our results are not signifi- 
cantly biased by this incompleteness at low cluster mass, 
we also show the cluster specific SFR versus M 2 oo for a 
sub-sample of clusters at z < 0.06 in the right panel of 
Figure [5] For these 22 low-z clusters, we are sensitive to 
star-forming galaxies down to SFR>1.4 MqJTT , which 
corresponds to the W4 band 5er detection limit of 6 mJy 
at z = 0.06, and is illustrated with a dotted curve. 

Using the sample of 62 clusters, the left panel of Fig- 
ure [5] shows no significant correlation between cluster 
specific SFR and cluster mass, over more than one or- 
der of magnitude in M2oo- The first pass coadd process- 
ing used here has the same calibration as in the WISE 
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Preliminary Data release, but was run every other day 
using only one day of data, creating gaps in coverage 
at low ecliptic latitude. These gaps should not signifi- 
cantly impact this result, since location of the gaps are 
random with respect to the cluster centers. We con- 
firm this expectation by examining the cluster specific 
SFRs of a sub-sample of 24 clusters with ecliptic lati- 
tude b > 44 deg. These clusters are at sufficiently high 
ecliptic latitude such that they have fairly complete spa- 
tial coverage with the first-pass WISE coadd data. The 
average cluster specific SFRs of the high ecliptic latitude 
sub-sample are overplottcd in Figure [5] in small purple 
diamonds. These high ecliptic latitude clusters show no 
significant correlation between cluster specific SFR and 
cluster mass, similar to the lack of correlation seen with 
the full cluster sample. 

The lack of correlation between cSSFR and M200 is 
also confirmed with the z < 0.06 sub-sample, for which 
we have applied a significantly lower detection limit of 
SFR>1.4 Moyr" 1 (Figure right panel). With the 
lower SFR limit, the clusters in Figure [5] have an average 
of 5.2 star-forming galaxies within R<3R200j with a large 
scatter. Five of the low redshift clusters have between 10 
to 16 star-forming galaxies. In contrast, for the sample of 
62 clusters with the demi-LIRG cut (SFR>4.7 Meyr" 1 ), 
there is an average of 1.6 star-forming galaxies per clus- 
ter within R<3R2oo- While the average number of star- 
forming galaxies per cluster is small in this case, the 
cSSFRs calculated using only the demi-LIRGs are well 
correlated with the cSSFRs using the SFR>1.4 M yr _1 
cut, with a linear Pearson correlation coefficient of 0.9. 
Therefore, even with an average of only 1.6 star- forming 
galaxies per cluster, the cSSFRs of the full cluster sample 
are qualitatively representative of cSSFRs that sample 
the fainter end of the infrared luminosity function. 

The binned data in both panels of Figure [5] show cluster 
specific SFRs that are consistent within la of each other, 
across the full range of M 2 oo ■ The error bars associated 
with the average cSSFR are the Poisson uncertainty on 
the number of star-forming galaxies detected above the 
SFR limit in each M200 bin. The lack of a significant 
trend between cSSFR and M200 indicates that trans- 
formation mechanisms which scale strongly with cluster 
mass may not play a dominant role in the evolution of 
star formation in clusters. 

While Figure [5] illustrates the integrated cluster SFR 
per M200 by taking the sum of all star-forming galaxies 
within R<3R2oo> our results remain robust when consid- 
ering only those galaxies within R<R2oo- Since galaxies 
within R200 are more susceptible to global cluster pro- 
cesses such as ram pressure, the lack of correlation be- 
tween cSSFR (R<R2oo) and cluster mass re-affirms the 
main conclusion drawn from Figured]- mechanisms that 
are strongly dependent on cluster mass do not play a sig- 
nificant role in the evolution of star formation in clusters. 

Two cluster processes that strongly scale with clus- 
ter mass are galaxy harassment and ram pressure. Ram 
pressure is a mechanism that acts to compress and/or 
strip cold gas reservoirs in galaxies. The impact of ram 
pressure is proportional to pv 2 , where p is the density 
of the intracluster medium (ICM) and v is the veloc- 
ity o f a galaxy with respect to the ICM ([Gunn fc Gotti 
1972). Galaxy harassment refers to the cumulative effect 
of high velocity close encounters between galaxies and 



has the potent ial to disturb mo rphologies and quench 
star formation ()Moore et al.lll996l) . The effects of galaxy 
harassment, like ram pressure, are expected to scale up 
with cluster mass, since both cluster velocity dispersion 
and ICM gas density increase with cluster mass. 

While there has been evidence in support of 
both ram pressure and harassment having an im- 
pact on star formation in individual cluster galax- 
ies fe.g. iRoediger fc Henslen [20051: iVollmer et aT]|200l 
iHaynes et al.l 20071 ) . the lack of a correlation between 
cluster specific SFR and cluster mass in Figure [5] indi- 
cates that neither of these mechanisms significantly trig- 
ger/quench star formation in local galaxy clusters, within 
the mass range of 10 14 to 7 x 10 14 M Q . 

O ur results are consiste nt with those of iGotol ()2005l ) 
and iPopesso et al.1 (|2007l ). who use a sample of ~100 
clusters at z < 0.1 and find no significant correlation be- 
tween Ha-derived clust er specific SFR and cl uster mass. 
Both IGotol (pOOTih and IPonesso et all (12001 use SDSS 
DR2 data, with the same limiting r-band magnitude as 
used in t his paper. The limiting H a SFR is ^3 M Q yr _1 . 
Similarly, IBalogh & McGed (|2010D demonstrate that the 
passive galaxy fraction is roughly constant as a function 
of system mass, from M ~ 10 13 M Q to 10 15 M . Assum- 
ing that the non-passive galaxy fraction is dominated by 
star-forming galaxies ra ther than AGN, the results of 
IBalogh &: McGeel (|2010l ) support our finding that cluster 
specific SFR is not strongly dependent on cluster mass. 

In contrast, several other studies have found evi- 
dence for an anti-c orrelation b e tween clust er specific SFR 
and c luster mass. iFinn et all ([20051 ) and iKovama et all 
(2010) obtain cluster SFRs by integrating Ha derived 
SFRs for members out to O.5R200, and cluster masses 
from observed velocity dispersions for a sample of eight 
and nine clusters, respectively. Their clusters span an 
approximate redshift range of 0.2 to 0.8, and the Ha 
SFRs are sensitive to ^1 M yr _1 . The anti-correlation 
observed between cluster specific SFR and cluster mass 
is significantly weakened when only clusters of the same 
epoch are considered, indicating that evolutionary effects 
may be the dominating factor. 

In addition to interactions between galaxies and the 
cluster environment, studies have demonstrated that 
tidal interactions between galaxies or galaxy-galaxy 
mergers can trigger a burst of star form ation (e.g. 
lOwen et alj [20051 iMartig fc Bournaudl [20081 ) . Galaxy- 
galaxy interactions are generally optimized in intermedi- 
ate density regions with low velocity dispersions, such as 
in galaxy groups, whereas the more massive clusters are 
less likely to host galaxy mergers. 

Figure [5] includes a substantial number of low mass sys- 
tems (57 clusters with M2oo< 3 x 10 14 M Q ), yet shows 
no sign of significantly elevated cluster specific SFR in 
low mass clusters relative to intermediate mass clusters. 
Any signature of enhanced star formation in low mass 
clusters due to galaxy-galaxy interactions should be par- 
ticularly evident in our data, because we integrate the 
SFR of cluster members out to three times the virial ra- 
dius, where galaxy-galaxy mergers and tidal interactions 
are more like to occur. A comparison of the distance to 
nearest neighbor among the demi-LIRG population and 
the general cluster population as a function of projected 
radius shows that the two populations have mean nearest 
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Fig. 5. — Left: Cluster specific SFR versus cluster mass for 62 clusters, with the average value per mass bin overplotted with large green 
square symbols. The average cluster specific SFR calculated using only clusters at high ecliptic latitude (and therefore greater coverage 
with WISE) are shown as purple diamonds, offset from the green square symbols for clarity. The subsample of clusters at high ecliptic 
latitude shows a similar trend relative to the full sample of clusters, indicating that the gaps present in the WISE coadd data at low 
ecliptic latitudes do not have a significant impact on our result. The dotted line shows the cluster specific SFR when SFR=4.6 M0yr — 1 
corresponding to the 6 mjy W4-band detection threshold at z = 0.1. Right: Cluster specific SFR versus cluster mass for a low redshift 
subsample of 22 clusters at z < 0.06. The dotted line shows a detection limit of SFR=1.4 MQyr - 1 , corresponding to the 6 mjy W4-band 
detection threshold at z = 0.06. There is no apparent trend between cluster specific SFR and M200 even with a lowered SFR threshold. 



neighbor distances that are consistent within la. This 
indicates that galaxy-galaxy mergers or close tidal inter- 
actions do not play a dominant role in the star formation 
properties of these clusters. 

Our results imply that even in the low mass ~ 5 x 10 13 
Mq groups, star formation has already been quenched to 
similar levels observed in clusters that are more massive 
by over an order of magnitude. The lack of correlation 
between cSSFR and M 2 oo (Figure [5]) , combined with the 
suppressed demi-LIRG fraction at 3R200 relative to the 
field (Figure01 would suggest a scenario in which star for- 
mation is quenched in a significant population of galaxies 
within small groups or filaments. Then once the remain- 
ing star-forming galaxies are accreted into the cluster en- 
vironment, a mechanism that is not strongly dependent 
on cluster mass must operate. We suggest strangulation 
as a plausible candidate for quenching star formation in 
cluster star-forming galaxies. 

As a galaxy enters the cluster ICM, its hot halo gas is 
removed, thereby cutting off its resource for future cold 
gas supplies, since the halo gas wo uld otherwise even- 
tually cool and se ttle onto the disk (Balogh et alJl2000t 
iBekki et al.l 120021) . The effectiveness of strangulation is 
less dependent on cluster m ass than ram pressure or ha - 
rassment. Simulations from lKawata fc Mulchaevl (2008) 
have shown that strangulation can occur in low mass 
groups, where the ICM-galaxy halo interaction is weak 
relative to clusters. Strangulation has been suggested by 
seve ral cluster studies as an important cluster me chanism 
(e.g. lMoran eTalllp^ |20^ iChung et al.ll2010D . includ- 
ing bv lTreu et al.l ( 2*0*0*31) who discovered "mild gradients" 
in the morphological fractions of galaxies in C10024+16 
out to large cluster radii, which were best explained by 
a slow-working gentle mechanism such as strangulation. 

While strangulation may quench star formation in 
massive clusters, our results suggest that the dominant 
process for regulating star formation in dense environ- 
ments does not depend upon cluster mass, and hence 



must be efficient in small groups or filaments. 

5. CONCLUSIONS 

We use data from WISE and SDSS DR7 to study 
the dependence of star formation on cluster mass and 
location within a cluster in 69 local (z < 0.1) galaxy 
clusters. Cluster membership is determined from SDSS 
DR7 spectroscopic redshifts for galaxies brighter than 
M r = —20.3, and star formation rates are determined 
with 22/im photometry from WISE using the relations 
outlined in lRieke et ail (120091 ) . 

Out of 69 clusters, we hnd a total of 109 star-forming 
demi-LIRGs with SFR> 4.6 M yr- X within 3R 200 - Both 
the fraction of demi-LIRGs and the mean specific SFR 
of cluster galaxies increases with projected distance from 
the cluster center. However, the fraction of demi-LIRGs 
remains below the field value even at 3 times the clus- 
ter virial radius. One plausible explanation for the sup- 
pressed demi-LIRG fraction at 3R200 is that even galax- 
ies that reside significantly beyond the virial radius have 
already been quenched in their previous environments, 
such as small groups or filaments. 

We also investigate the impact of cluster mass on star 
formation by presenting the total cluster SFR normalized 
by M200 as a function of M2oo- We find no evidence of a 
correlation between the cluster specific SFR and cluster 
mass in this first uniform dataset to detect obscured star 
formation out to several times the virial radius in a large 
sample of low redshift clusters. Our result indicates that 
cluster mechanisms which scale with cluster mass, such 
as ram pressure or harassment, are not likely to play a 
dominant role in the evolution of star formation in local 
clusters. 
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TABLE 1 List of demi-LIRGs from 69 clusters sorted in descending order of 
L IH . Galaxies marked as AGN have WISE color Wl-W2>0.5 and are 
assumed to have a mid-IR flux dominated by an AGN rather than star 
formation. The six BPT Seyferts with W1-W2 color consistent with star 
formation rather than AGN, are marked with an asterisk. 
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